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Abstract 

Using 420 pb" 1 of data collected on the T(55) resonance with the CLEO III detector, we recon- 
struct B mesons in 25 exclusive decay channels to measure or set upper limits on the decay rate 
of T(5S I ) into B meson final states. We measure the inclusive B cross-section to be a(T(5S) — > 
BB(X)) = (0.177 ± 0.030 ± 0.016) nb and make the first measurements of the production rates of 
<t(T(55) -► B*B*) = (0.131 ± 0.025 ± 0.014) nb and cr(T(5S") -»• BB*) = (0.043 ± 0.016 ± 0.006) 
nb, respectively. We set 90% confidence level limits of ct(T(5S) -► BB) < 0.038 nb, <r(T(5S) -> 
B(*)b(*)tt) < 0.055 nb and a(T(5S) -> BBtttt) < 0.024 nb. We also extract the most precise value 
of the B* mass to date, M(B*) = (5411.7 ± 1.6 ± 0.6) MeV/c 2 . 
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The T(55) resonance was discovered by the CLEO [0] and CUSB 0] collabora- 
tions. Its production cross-section and mass were measured to be about 0.35 nb and 
(10.865±0.008) GeV/c 2 Q ^respectively. Final states can be: BB, BB*, B*B*, BBix, BB*n, 
B*B*n, BBixix, B S B S , B S B* and B*B*. Here, B=B U or B d , and the 7r may be charged or 
neutral (consistent with charge zero of the final state). Throughout this article, we use BB* 
to signify both BB* and B*B. Including a symbol in parentheses indicates that it may or 
may not be present. The B cross-section in this region is well-described by the Unitarized 
Quark Model (UQM) 0, which predicts that about 1/3 of the bb decay rate is to B^B^ 
and that B*B* dominates the inclusive B rate. A previous CLEO measurement using inclu- 
sive D s production revealed that T(5S) -> B^B^ constitutes (16.0±2.6±5.8)% of the total 
bb rate j^]. A second analysis which performed exclusive reconstruction of B s mesons 
found a(e + e~ -> B*B* S ) = (0.11±°:o3 ± 0.02) nb (about 1/3 of the total hadronic resonance 
cross-section). The two results are consistent with each other and with predictions of the 
UQM. 

In this Letter, we measure the contributions of various B meson final states to the T(5S) 
decay. These measurements may better constrain coupled-channel models in the T mass 
region as well as near the lower ip resonances p. We also exploit exclusively reconstructed 
B mesons from this analysis and the corresponding B s analysis [5] to extract the most precise 
measurement of the B* mass to date. 

CLEO III is a general purpose solenoidal detector that includes a tracking system for 
measuring momenta and specific ionization (dE/dx) of charged particles, a Ring Imaging 
Cherenkov detector (RICH) to aid in particle identification, a Csl calorimeter for detection 
of electromagnetic showers, and a muon system for identifying muons 0]. 

The analysis presented here uses 420 pb -1 of data collected on the T(5S) resonance 
(y/s = 10.868 GeV) at the Cornell Electron Storage Ring. Using techniques pioneered at 
the T(AS), we utilize two kinematic variables: the energy difference AE = E^ eSum — Eb and 



the beam-constrained mass Mb c = y E^ egm — p\, where E B (p B ) is the energy (momentum) 
of the reconstructed B candidate and -Ebeam is the beam energy. Because of its low energy, 
reconstruction of the photon in B* — > B^ is not essential; to maintain high efficiency, we do 
not reconstruct the B* . (Charge conjugate final states are implied throughout this Letter.) 

To obtain a B meson sample of high purity, events are required to have at least five 
charged tracks and a ratio of the second to zeroth Fox- Wolfram moment Q , i?2 < 0.25. 
Candidate B mesons are reconstructed in exclusive final states containing either J/ip or 
mesons. 

Charged particles are required to pass standard selection criteria and are identified by 
using their measured momenta in conjunction with dE/dx, RICH, calorimeter and muon 
system information. For particle types i, j (i, j = n,K,p) we define x 2 -like quantities for 
dE/dx as the difference in the measured and expected dE/dx, normalized by the uncertainty, 
i.e., xf /dx = (dE/dxf 1 ™ 8 - dE/dxf v )/o % , and for RICH as Ax?j = d - Cj (difference 
in negative log-likelihood between hypotheses i and j), respectively. We require at least 3 
detected Cherenkov photons from the RICH. Pions are identified by requiring \xt E ^ dx \ < 4 or 

AxI^k < 5. For kaons, we define a combined quantity, Xcomb = {Xk^^ 2 ~ {xt E ^ dx ) 2 + ^Xk^ 
and require Xcomb < 0- Electron candidates are formed from particles that have a ratio of 
calorimeter energy (E e ) to measured momentum (p e ) in the range 0.5 < E e /p e < 1.25. 
Muons are identified by either having penetrated at least 3 layers of iron absorber or by 
having deposited energy in the calorimeter consistent with a minimum ionizing particle 
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(E < 300 MeV). Photons are formed from showers that have deposited at least 30 MeV of 
energy in the calorimeter and are not associated with a charged track. Pairs of photons that 
have an invariant mass within 2 standard deviations (a ~ 6 MeV/c 2 ) of the known 7r° mass 
(Mn-o) [^| are defined as 7r° candidates and are kinematically constrained to give M n o. 

Candidate J/^'s are formed from p + p~ or e + e~ pairs. For muon pairs, we require 
3.05 < M M+M - < 3.14 GeV/c 2 . For e+e" combinations with 1.50 < M e+e - < 3.14 GeV/c 2 , 
bremsstrahlung photons are searched for among the showers with no matching charged track 
and within a 5° cone about each electron's initial direction. For each p + p~ and e + (7)e~(7) 
candidate, we perform a mass-constrained fit to the J/ip mass P and make a loose re- 
quirement that the fit x 2 P er degree of freedom is less than 100. Candidate p + (K s ) [K*°] 
mesons are formed from tt + 7t° (7r + 7r~) [f^ + 7r~] combinations that have an invariant mass 
in the range from 620-920 (490-505) [820-970] MeV/c 2 . D + (D°) meson candidates are 
reconstructed via their decays to K~ir + 7i + (K~tt + , K^7i + tt + ti^ and K~7t + 7t ) and are re- 
quired to have an invariant mass within 2cr of their PDG !9j values. To reduce combinatorial 
background in D° — ► K^tt+tt , we require p n o > 400 MeV/c. Candidate D* + — > D°7r + 
(D* + — > D + n°) decays are formed from D and tt candidates that have a mass difference 
in the range 144 < M D *+ - M D o < 147 MeV/c 2 (139 < M D *+ - M D+ < 143 MeV/c 2 ). 
Similarly, D*° mesons are reconstructed in D 0/ k°, and the mass difference is required to be 
in the interval 140 < M D *o - M D o < 144 MeV/c 2 . 

Candidate B mesons are reconstructed in the 25 decay channels listed in Table |U For 
B -> Dp and B -> L>*p 0, we take advantage of the helicity angle (6^) ^l| distribution 
in these decays and require | cos 9h\ > 0.3. To improve the signal-to-background ratio, we 
also reject low momentum (backward-emitted) 7r°'s from the p + decays by requiring cos 8h > 
—0.7. Table U also shows the product branching fractions, £>j, including the branching ratios 
of the daughter modes [9|, and the reconstruction efficiencies, determined from Monte Carlo 
simulations 0, 0, Q] of these decays followed by a geant 15] based detector simulation. 
We validate our simulation and analysis procedure by measuring branching fractions for 
these decay modes using data collected on the T(4S) resonance. Good agreement with the 
world averages are found for all modes. 

We first determine the total B meson yield by fitting the invariant mass distribution 
formed from candidates in the M bc , AE region of 5.272 < M bc < 5.448 GeV/c 2 , -0.2 < 
AE < 0.45 GeV. The relatively wide AE region is used to avoid biasing the background 
shape. The invariant mass distribution, shown in Fig. Q is fit to the sum of a second-order 
polynomial background and a Gaussian signal shape whose width is fixed to 12.3 MeV/c 2 , 
the expected average resolution of these candidates. We find a yield of 53.2 ± 9.0 events; 
fitting to the J/ip and Z?W distributions individually results in 11.2±3.5 and 42.3±8.4 
events, respectively. Using J2&i e i — 7.2 x 10~ 4 (see Table |TJ), we measure a cross-section 
a(T(5S) -> BB(X)) = (0.177 ± 0.030) nb. 

Figure El shows the reconstructed events in the M^ c — AE plane for T(5S) data. Signal 
and sideband regions are defined using MC simulations of these final states. To extract rates 
for BB, BB*, and B*B* separately, we select events in a signal region defined by the area 
between the diagonals M bc = 1.018 AS + 5.248 GeV/c 2 and M bc = 1.018AE + 5.312 GeV/c 2 . 
This restricted signal region has a total Y^^i e i — 5.7 x 10~ 4 . Lower and upper sidebands 
of the same AE width, also shown in Fig. |21 are shifted to the left and right of the signal 
region by 10 MeV, respectively. 

The BB, BB* and B*B* final states are kinematically well separated, but B^B^n final 
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TABLE I: Modes used in exclusive B meson reconstruction along with their product branching 
fractions (£>j) Q, reconstruction efficiencies (e^), and expected yields (iV* ) in data assuming 
a(T(5S) -» BB(X)) = 0.2 nb. 
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states have a large degree of overlap, and with the limited statistics cannot be distinguished. 
The BBnn final states, because of the limited phase space, peak at M^ c ~ -Ebeam- If their 
rate is large enough, their shape (in M^ c ) will be sufficient to distinguish them from the 
broad tail of B^*'B^*'tt final states that extend into the M^ c region of BBttit. 

Events in the signal region of Fig. |21 are projected onto the M^ c axis (see Fig. |3J) and 
fit to the sum of a flat background and three Gaussian signals, one each for the M^ c peaks 
produced by BB, BB* and B*B* events. The signal resolutions are set to o =4.0, 6.2 and 
7.0 MeV/c 2 , respectively, as determined from MC simulation, and the background is fixed 
to 0.7 events/4 MeV, as determined from the average of the upper and lower sidebands. In 
the fit we use the precisely known mass difference Mb* — Mb j9| and constrain its value to 
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FIG. 1: The B meson invariant mass distribution for all B decay modes listed in Table [I] in T(5S) 
data. The points are the data and the curve is the fit described in the text. 

47.5 MeV/c 2 ^6|. The middle peak, which corresponds to T(55) — > BB*, is not constrained 
in the fit and is found to be within la of the expected value. 

The fitted yields are 3.7±|i BB, 10.3 ± 3.9 BB* and 31.4±6.1 B*B* events. Only the 
latter two are statistically significant with significances, determined from the change in log- 
likelihood when the contribution from each peak is removed, of 4.3<r and 7.6a, respectively. 
For BB, we compute an upper limit of 7.5 events at 90% confidence level (CL). A potential 
excess m #(*)£(*) 7i is examined by plotting the invariant mass of candidates in the B^B^tt 
region defined by 5.351 < M bc < 5.429 GeV/c 2 and -0.2 < AE < 0.45 GeV, which should 
exhibit a peak at Mb (see inset in Fig. EJ). This Mb c — AE region includes (88 ±6)% of 
reconstructed B^*'B^*'tt events, where the uncertainty reflects the maximum variation based 
on the possible B^*'B^*>tt final states. That distribution is fit to the sum of a Gaussian signal 
whose mean and r.m.s width are constrained to 5.279 GeV/c 2 and 12.3 MeV/c 2 , respectively, 
and a linear background shape. The yield of 6. 7t^_l is not statistically significant, and we 
compute an upper limit of 13.1 events at 90% CL. 

For the BBtttt final state, we select events in the region 5.429 < M^ c < E^ eam GeV/c 2 
and —0.2 < AE < 0.45 GeV and find three events consistent with Mg. This additional 
requirement on M^ c and AE has an efficiency of (95 ± 2)%. While the combinatorial back- 
ground is small (~ 0.3 events), the cross-feed from B^*'B^tt into the BBtttt signal region is 
(12 ± 6)%. If we conservatively assume a B^'B^'tt yield equal to its 90% upper limit value 
and that the 3 BBtttt candidates are also B^B^tt, we would expect 1.0 - 2.9 B^B^tt 
events to lie within the BBtttt mass region. Based on this range of expected background 
and 3 observed events, we take the most conservative upper limit on T(5S) — > BBtttt, which 
corresponds to 6.4 17| events at 90% CL. For illustration, we superimpose on Fig. 6.7 
B^B^tt (lightly shaded, with a ratio of 1:1:1 ratio for B*B*tt:BB*tt:BBtt) and 3 BBtttt 
(darker shading) events. Yields, efficiencies, cross-sections and relative production fractions 
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FIG. 2: Scatter plots of M^ c vs. AE for all B decay modes listed in Table [I] in T(55) data. The 
diagonal lines show the expected signal (solid) and 2 sideband (dashed) regions, as discussed in the 
text. The horizontal lines show the regions for the various BB(X) final states. 



1630106-002 




FIG. 3: Distribution of M^ c for all reconstructed B modes listed in Tabled for T(5S) data. The 
histogram displays the data, the curve shows the fit described in the text and the flat line shows 
the background as determined from a fit to the sidebands as discussed in the text. Distributions for 
B^B^tt (lightly shaded) and BBtttt (darker shading) obtained from MC simulation of these final 
states are superimposed for illustrative purposes. The inset shows the invariant mass of candidates 
m the B^B^n M bc region with the fit superimposed. 
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are summarized in Table |HI We also show the cross- sect ions as determined from the J / ip 
and modes separately. We find that B*B* is indeed dominant, comprising (74 ± 15)% 
of the BB{X) rate. 

Several sources of systematic uncertainty on the cross- sect ions measurements are consid- 
ered. Potential errors from the background normalization and shape are evaluated by using 
different background parameterizations and varying the normalization within its uncertainty. 
The corresponding uncertainties in the cross-sections vary from 3.1% for T(5S*) — > BB(X) 
to 16.7% for BB. Uncertainties in the reconstruction efficiencies include contributions from 
charged particle tracking and identification, and 7r° reconstruction, and finite MC statis- 
tics. Averaged over all modes, we find an uncertainty of 6.5%. The analysis procedure was 
also checked by comparing I?-meson branching fractions in our signal modes measured using 
data collected on the T(45) resonance with PDG values fsj. We find a relative difference of 
(1±3)%, averaged over all modes, indicating that the efficiencies are well understood. Errors 
due to the fixed signal shape parameters are determined by varying them within their un- 
certainties and refitting (3%-4%). The occurrence of multiple candidates in data (in a single 
event) are found to agree with simulation to within 3%. Uncertainties on input branching 
fractions and measured integrated luminosity contribute 3% and 2%, respectively. These 
systematic uncertainties are added in quadrature and the resulting values are included in 
the cross-sections shown in Table |H] 



TABLE II: Summary of yields, efficiencies, cross-sections and fractional contributions of various 
subprocesses in a(T(5S) — > BB(X)) decays. Upper limits are set at the 90% CL. Uncertainties 
are from statistical and systematic sources, respectively. 
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< 0.055 
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6.8 ±0.5 


< 0.024 
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a(T(5S) -> BB{X)) 


53.2±9.1 


7.2 ±0.5 


0.177 ±0.030 ±0.016 




J/tp Modes 


11.2±3.5 


6.3 ±0.5 


0.295 ± 0.092 ± 0.028 




flW Modes 


42.3±8.4 


0.9 ±0.1 


0.161 ±0.032 ±0.015 





We proceed to use the M^ c distribution from this analysis in combination with the one 
for B* in Ref. (if to obtain an improved measurement of the B* mass. Since those results 
used exactly the same data set as in this analysis, the largest systematic error, the beam 
energy calibration of (±4.6 ± 2.9) MeV p, cancels out in the (uncorrected) M^ c difference, 
M hc (B*) - M hc (B*). The rightmost peak in Fig. 01 corresponds to T(55) -> B*B*, and 
its mean value is measured to be (5333.1 ± 1.3(stat)) MeV/c 2 . The M bc peak value for 
T(55) -> B*B* from Ref. 18] is (5418.2 ± 1.0 ± 3.0) MeV, where we have added back the 
(±4.6 ± 2.9) MeV beam energy correction to obtain an uncorrected value. The difference in 
the Mb c peak values, M^, C (B*B*) — M\, C (B*B*), can be translated into the mass difference, 
M(B* S ) - M(B*) after correcting for the -1.7 (-0.1) MeV/c 2 bias that is introduced due to 
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our use of reconstructed instead of Bf-, mesons. We therefore find a mass difference 
M(B* S ) - M(B*) = (86.7 ± 1.6 ± 0.2) MeV/c 2 . The 1.6 MeV/c 2 error is statistical and the 
0.2 MeV/c 2 uncertainty is from systematic errors in fitting pur M^ c spectrum. Combining 
this mass difference with M(B*) = (5325.0 ± 0.6) MeV/c 2 |9j, we obtain an improved value 
for the B* mass, M(B* S ) = (5411.7 ± 1.6 ± 0.6) MeV/c 2 . Using the well-measured B s mass 
from CDF of M(B S ) = (5366.01 ± 0.73 ± 0.33) MeV/c 2 Q, we determine the 1~ - (T mass 
splitting M(B*) - M{B S ) = (45.7 ± 1.7 ± 0.7) MeV/c 2 . This mass splitting measurement 
supersedes the previous CLEO result of (48 ± 1 ± 3) MeV/c 2 @ and is significantly more 
precise than an earlier indirect measurement of (47.0 ±2.6) MeV/c 2 It is also consistent 
with the corresponding splitting in the Bd system of (45.78±0.35) MeV/c 2 f9j as expected 



from heavy-quark symmetry [2l| 



In summary, we have measured or set upper limits on the rates for the various final states 
in T(5S) decay. We find that predictions of the UQM j3[ are consistent with our findings 
that B*B* is dominant, with a measured value of (74 ± 15)% of the total B rate. The 
BB* rate is measured to be about 1/3 of the B*B* rate. Upper limits on BB, B^B^tt 
and BBtttt have also been presented. Lastly, we utilized the M^ c peak positions for B* and 
B* to extract M(B*) = (5411.7 ± 1.6 ± 0.6) MeV/c 2 , which is the most precise value of 
the B* mass to date. 
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